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The color-suppressed D°uj decay modes have just been observed for the 

first time. The rates are all larger than expected, hinting at the presence of final state interactions. 
Considering ^ mode alone, an elastic D^*\ —> dm TT rescattering phase difference 

S = Si /2 — S 3/2 ~ 30° would suffice, but the D^*^°r], D^uj modes compel one to extend the 

elastic formalism to SU(3) symmetry. We find that a universal a^jax = 0.25 and two strong phase 
differences 20° ~ 0 < 5 < 5' ~ 50° can describe both DP and D*P modes rather well; the large 
phase of order 50° is needed to account for the strength of both the and D^*^° r] modes. For 

DV modes, the nonet symmetry reduces the number of physical phases to just one, giving better 
predictive power. Two solutions are found. We predict the rates of the ^ Df K~, D*'^K~, 
D°p^, Df K*~ and modes, as well as ^ D°K^, D*°K°, D°K*° modes. The formalism 
may have implications for rates and CP asymmetries of charmless modes. 

PACS numbers: 11.30.Hv, 13.25.Hw, 14.40.Nd 


I. INTRODUCTION 

The Belle Collaboration has recently observed [Q the 
D^rj and decay modes, as well 
as finding evidence for 5° ^ D*'^ri and D*°uj. The decay 
branching ratios {B) are all at a few times 10“"^ level: 
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The CLEO Collaboration has also reported the obser¬ 
vation of —!■ modes, 

B{W D°Tr°) = (2.74 toU ± 0.55) x 10-^ 
B{W D*°Tr°) = (2.20 ± 0.79) x 10-^ 

with rates in agreement with Belle. These modes are 
usually called color-suppressed B decays. In contrast 
to the much faster “color-allowed” 5° ^ D'^7r~ {B ~ 
3 X 10“^ [Q) decay where 7r“ is emitted by the charged 
current (Fig. 1(a)), there is a color mismatch in forming 
the meson from c and u produced by & ^ cud decay 
(Fig. 1(b)). We have indicated in Figs. 1(a) and (b) the 
effective Wilson coefficients, ai and 02 [||, that is respon¬ 
sible for the decay. For the charged B~ —> D^-k~ decay, 
both ai and 02 type of diagrams contribute. 

The factorization of the B^ —s- decay amplitude 

has recently been demonstrated to follow from QCD in 
the heavy b quark limit & n|, and oi as computed from 
QCD factorization is close to 1-1.1 from “gener¬ 

alized factorization” Q. On the other hand, it is known 


that the effective coefficient cannot be calculated 
in QCD factorization ||^. It is remarkable that the 
value as extracted from the celebrated B 
modes agrees rather well with the value extracted from 
B~ D^tt~, i.e. consistent with ~ 0.2-0.3 and 
real. However, the rates as observed by Belle 

and CLEO are considerably higher than the general¬ 
ized factorization estimates Q using this value, 

suggesting that final state rescatterings (FSI) such as 
£)(*)+ 7 r“ ^ could be active. Alternatively, it 

could indicate non-universal 02 , that is, 02 for color- 
suppressed modes are larger than i 02 for J/ijjK^*^ 
modes, and are furthermore complex to accommodate 
B- TT . Whether in the form of FSI rescattering 

or a larger and complex 02 , one in general acquires com¬ 
plexity of decay amplitudes in the form of CP conserving 
phases. This may have implications for direct CP viola¬ 
tion asymmetries, since similar effects may be present in 
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FIG. 1: Color (a) allowed and (b) suppressed B Dtt decays. 
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processes such as charmless B Ktt modes. 

The experimental data in fact compels one to broaden 
the horizon. The Belle results on D^rj and 

are considerably higher than generalized factorization es¬ 
timates 0 as well, suggesting that one needs to go be¬ 
yond Dtt —> Dtt considerations in the FSI framework. 
The alternative of having process dependent a 2 ’s ||^, ||] 
would imply loss of predictive power. While it is clear 
that 02 in general will be process dependent, we do ab¬ 
hor the loss of predictive power. In this paper we wish 
to explore how the situation could be remedied. 

Since the data is new and still rather incomplete, our 
approach will be phenomenological, without aim for rigor 
or completeness. Let us start from the isospin decompo¬ 
sition of B —> D^7r~, D~^7r~ and decay amplitudes 
that is lucidly outlined in Ref. Q, 


Adott- — ■v/3^3/2) 



where the final state is emphasized. In the absence of 
FSI, the smallness of Z)° 7 r° rate can be viewed as due 
to cancellation between the Axj^ and ^ 3/2 amplitudes, 
which are real under factorization. But these amplitudes 
in general become complex under FSI and Aj^Ot^o stands 
to gain strength. However, the isospin or triangular re¬ 
lation Ajj+t^- = ^/2Ajjo^o + Ajjot^- always holds. We 
note that Ajjot^- is purely H. 3/2 and can only rescat¬ 
ter into itself. It is therefore reasonable to maintain 
as extracted traditionally, i.e. the same as from 
since it is too good to be just a coinci¬ 
dence. 

The formalism of Eq. (||) can generate B^ 
by elastic rescattering from the color-allowed —> 

TT mode [IL The problem is the strength of the 
B^ —> and ^ modes observed by Belle. 

The ?7 and uj are isosinglets, and the 77 mass is quite dif¬ 
ferent from pions. Thus, to generate and 

final states from and p~, strictly speak¬ 

ing, involves inelastic rescattering, although only in the 
isospin 1/2 channel. However, it may be reasonable to 
extend from isospin (an excellent symmetry of the strong 
interaction) to SU(3), expecting that the latter becomes 
a good symmetry at the ttib rescattering scale. 

Let us qualify the last statement further. SU(3) fla¬ 
vor symmetry is a symmetry of the strong interaction 
as QCD is flavor blind, except that the flavor symme¬ 
try is broken by quark masses. Thus, in terms of SU(3) 
multiplets, masses vary within the multiplet according to 
SU(3) breaking, and meson production differs in strength 
as reflected in, for example, the decay constants and tran¬ 
sition form factors. But if we enlarge the isospin doublet 
£)(*) = (£)(*)+, £)(*)0) and triplet tt = (tt^, tt®, 7r“) 
to SU(3) triplet = (£»(*)+, and the 

meson octet H (likewise from p to vector V), we note 


that rescattering occurs at ms Wq scale, hence the 
DH ^ DH strong rescattering amplitude should respect 
SU(3) symmetry to good degree. SU(3) breaking effects 
are taken into account in the initial meson formations 
from B meson weak decay, which is done in the (QCD) 
factorization framework with and effective 

short distance Wilson coefficients. 

Thus, our aim is to extend the usual elastic Dn —> Dtt 
rescattering to guasi-elastic DH ^ DH rescattering in 
final state, so and io modes are naturally in¬ 

corporated. Our framework is rather close in spirit to 
the original isospin analysis of H ^ Dtt, and is as close 
to elastic as one can get. The most general formalism in¬ 
volving inelastic rescatterings from all possible hadronic 
final states cannot be tackled. In general it involves large 
cancellations and, statistically speaking, small phases 
Hopefully, and in a sense true by duality, the inelastic ef¬ 
fects are contained already in and a^***^. 

In Sec. H we introduce the general framework of a T- 
matrix, discuss its link to the optical theorem, and also 
fix the phase convention. The formalism is applied to 
the B Dtt modes in Sec. HI. Three types of rescatter¬ 
ing amplitudes are identified: a diagonal “Pomeron”-like 
piece, and two “Regge”-like pieces denoted as “charge- 
exchange” and “annihilation”. We relate these to the 
usual isospin 1/2 and 3/2 rescattering phases, and show 
that only 5 = 81/2 ~ <^ 3/2 matters, as expected. 

In Sec. IV we extend from SU(2) to SU(3) multi¬ 
plets in the final state. For DP and D*P modes, where 
P stands for the pseudoscalar octet, the extension is 
straightforward. The question of whether to include the 
flavor singlet rji in the final state is bypassed by noting 
a) the absence of data, which would remain the case for 
a while unless B'^ D^rj' turns out to be much larger 
than ^ D^p, and b) Ua( 1 ) anomaly that singles out 
the 771 field, and perhaps as a consequence, c) relatively 
small singlet-octet or 77 - 77 ' mixing, allowing us to identify 
778 ~ 77 for convenience. We thus ignore 771 completely in 
this work. We find the previous picture of three types of 
rescattering parameters still hold, but one now has two 
i5-like phase differences. The extension to DV final states 
is treated differently. By noting that the vector mesons 
satisfy U(3) rather than SU(3) symmetry, we use a nonet 
V held rather than an octet one. We refrain from dis¬ 
cussing D*V modes since data is scarce, and since two 
helicity (or partial wave) amplitudes are involved. 

We carry out a numerical study in Sec. V. For the 
DP modes, we hnd two sets of solutions for the two FSI 
phase differences, which are of order 20° and 50°. One 
solution is similar to the D*P case, and has a very tiny 
B^ Df K~ decay due to the smallness of “annihi¬ 
lation” rescattering. The other solution gives B{B^ 
D+K-) ~ 5 X 10 which is ruled out by experiment. 
For DV modes, we have two solutions: one does not 
have annihilation contribution hence B{B^ D^p^) ~ 
H(H° ^ D^ui), while the other does not have exchange 
contribution and B{B'^ D°p°) > B{B^ D^lo^). 
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In Sec. VI we compare our plus FSI rescatter¬ 

ing” approach with the viewpoint of “process dependent 
02 ”, and discuss possible future applications. The con¬ 
clusion is then offered, followed by Appendices that give 
same results from a SU(3) decomposition approach and a 
geometric (triangular) representation of our rescattering 
picture. 



FIG. 2: Illustration of optical theorem, Eq. (^. 


II. FINAL STATE RESCATTERING 
FRAMEWORK 


Before we turn to B ^ Dir decay applications, let us 
give the isospin structure of the related amplitudes. The 
responsible effective weak Hamiltonian is given by 


Let iLw denote the weak decay Hamiltonian. We 
assume absence of weak phases (or they are factored 
out), hence, from time reversal invariance of 7Lw(= 
UtH^U^), one has, 

{i-ont\H^\BY =Y,S*S\ont\H^\B), (2) 
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Hw = -^I4bl4d[cifi7^(l-75)'UC7;,(l-75)6 

-hC 2 C 7 '"(l- 75 )'ud 7 ;,(l- 75 ) 6 ], (7) 

where VcbV*^ can be treated as real for our purpose. It 
is clear that Hw transforms like a / = 1 , J 2 = 1 vector 
under isospin Q. It is useful to state explicitly the phase 
convention and the isospin structure of these mesons: 


where Sij = (f;out|j;in) is the strong interaction S- 
matrix element, and we have used t/rlout (in))* = 

|in (out)) which also fixes the phase convention. Eq. (H) 
can be solved by (see, for example |^) 

(*;out|Hw|i?)=5]4^"A°, (3) 

i 

where is a real amplitude. To show that this is indeed 
a solution of Eq. (^, one needs to use = Sji, which 
follows from time reversal invariance of the strong inter¬ 
actions and the phase convention we have adopted. The 
weak decay amplitude picks up strong scattering phases 
Also note that since 5^/^ is unitary, we must have 

^|(*;out|Hw|i?)P=^|A°p. (4) 

i I 

Eq. (^ implies an identity related to the optical the¬ 
orem. Noting that S' = 1 -I- iT, we find 

2Im(z;out|Hw|S) = ^ r7(j; out|Hw|S). (5) 

3 

Thus, for B decay to a two body final state with momen¬ 
tum (pi,P 2 ), one has the relation 

2ImM(pB ^ P 1 P 2 ) 

y.M{pB ^ {qk})M*{pip 2 ^ {qk}), (6) 

which relates the imaginary part of the two body decay 

amplitude to the sum over all possible B decay final states 
{qk}, followed by {qk} P 1 P 2 rescattering. This equa¬ 
tion is consistent with the optical theorem to all orders 
of the strong interactions but only to first order of the 
weak decay vertex, as we illustrate in Fig. 2. 


\W) = 

hd) - 

bd\0) 

\B-) = 

hu) ^ 

Wo) 

\D+) = 

cd) ~ 

cd| 0 ) 

|H°) = 

cu) ~ 

CM 0) 
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-| 1 , 0 ). ( 8 ) 


Note that the isospin structure is defined according to the 
fields that create the states, such that they conform with 
the definition of the isospin structure of Hw- One could 
alternatively define isospin quantum numbers according 
to states and modify that of Hw accordingly. 

By isospin decomposition, we have Q (cf. Eq. (0)), 


Ad+t:- = ;out|Hw|.B°) 


+ y'^A3/2, 

(9) 

(H° 7 r°;out|Hw| 5 °) 

y|kIi/2 - y|7l3/2, 

( 10 ) 

(H°7r ;out Hw H ) 

V 5 A 3 / 2 , 

( 11 ) 


where 

Ai/2 = ((-D7r)i/2;out|Hw|H“) , 

^ 3/2 = ((-D7r)3/2;out|Hw|H°) 

= ■^((A’7r)3/2;out|Hw|H-). (12) 
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The last step follows from the Wigner-Eckart theorem. 
Eqs. (P)-(|rT|) imply the triangular isospin relation 


- V2Ad0j,.0 + . 


(13) 


We note that the sign of Ajjo^o in Eqs. (W) and M is 
different from Ref. Q , but is consistent with Ref. |^. It is 
easy to see from above equations that 1 ^ 1/2 P + 1 ^ 3/2 P = 

The isospin relations Eqs. (|^)-(|l^) are valid whether 
one has (in)elastic ESI or not. For example, assuming 
factorization hence ignoring ESI, one has 


Ado-itO ^DO^o — 

^DO-K- = T + C, 

Ao+n- ^ Af^^^_ = T + E, (14) 


where T, C, E are the color-allowed external lE- 
emission, color-suppressed internal lE-emission and IE- 
exchange amplitudes, which we shall discuss later. These 
amplitudes clearly satisfy Eq. (|l3). 

The general validity of Eqs. (|)-(^^ in fact allows 
one to extract Axii and ^ 3/2 directly from the measured 
TT rates without any further recourse to theory. Using 
the Belle and CLEO average of = (2.9 ± 

0.5) X 10-4 and B(^ = (2.5 ± 0.7) x IQ-^, 

with other rates and = 1.073±0.027 taken 

from PDG Q, we find 


' ^ ' - = 0.71 ±0.11 (0.75 ±0.08), 

y2|A3/2h’'P*- 

|^|expt. ^ 29° ±6° (30° ±7°), (15) 

for modes, where <5 = ^ 1/2 — < 13/2 is the phase 

difference between Ay ^2 and ^ 3 / 2 - This strongly suggests 
the presence of ESI Q, ||, . 

As we turn on ESI, although the isospin relations still 
hold, we would clearly lose control if the full structure 
shown in Eq. (^ is employed. Even if all possible B 
decay rates can be measured, it would be impossible to 
know the phases of each amplitude. Furthermore, we 
know very little about the strong rescattering amplitudes. 
However, the subset of two body final states that may 
be reached via elastic rescatterings stand out compared 
to inelastic channels. It has been shown from duality 
arguments [Q as well as a statistical approach Q that 
inelastic ESI amplitudes tend to cancel each other and 
lead to small ESI phases. We shall therefore separate 
{qk} into two body elastic channels plus the rest. We 
first explore the familiar Q H ^ D+tt-, D°'k~ and 
case, then try to stretch the scope of elasticity. 


III. ELASTIC ESI IN THE Dn SYSTEM 


factorization amplitudes of Eq. (0) , one has 


( AdOtt- 
Ad+tt- 
Ado-^o 


= 54/2 I A 

A 


A^ 

r-- 


r 


(16) 


DOttO 


A major assumption is involved here: we assume that 
one can separate hard and soft effects. The factoriza¬ 
tion amplitudes sum over “hard” contributions, includ¬ 
ing attempts at incorporating the effects of the largely 
intractable full set (generalized to n body) of inelastic 
amplitudes, illustrated for two body final state in Fig. 2. 
The matrix describes the nonperturbative ESI from 
the factorized Dtt “source” amplitudes. We illustrate this 
in Fig. 3, where the “hard” part is shrunk to a point, and 
we focus on two body elastic ESI. 

It is instructive to show how can be obtained. 

In the usual approach, one notes that with only elastic 
E>7r Dtt rescatterings, the S-matrix is diagonal in the 
isospin basis (isospin invariance), that is 

5diag = = diag(S'|_|, Si i), (17) 



Unitarity then implies that the diagonal elements of S 
can only be pure phases, or 

1 = 53 3 = (19) 

2 ’ 2 2 ’ 2 ^ 


54 /^ is likewise diagonal, i.e. (5^/^)i/2.i/2 = and 

( 54 /^) 3 / 2,372 = Elastic ESI is equivalent to A/ = 

g2(5/|^/| |^/| taken from factorization approach [|j. 

The isospin relation of Eq. (O) is clearly satisfied. 

For later use we express 5±^ in the basis of Eq. 0), 


/I 0 0 \ 

5i/^ = e*‘^3/2 0 i(l±2e*‘^) -^(l-e'-^) , 

Vo -^(1-e*^) i(2±e*^) / 

( 20 ) 

with overall phase (i.e. of Aua^^-) factored out, 

and 6 is the phase difference between isospin 1/2 and 
3/2 amplitudes, which is physical. 

An alternative way to obtain 5^/^ is through the opti¬ 
cal theorem, i.e. Eqs. (||) and (^). This approach is less 



Let us consider elastic final state rescattering in H ^ r-^ 

Dtt modes. By using Eq. 0, with A° taken as the lustration of Eq. (j^): two body rescattering. 
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familiar and more awkward than the previous one. How¬ 
ever, it has the advantage of being readily generalizable 
to encompass other modes such as 5° ^ D^rj, Df K~, 
which are now of interest and also provide physical in¬ 
terpretation. To use the optical theorem, we need to 
study the strong S matrix, or equivalently the T ma¬ 
trix. The T-operator corresponding to the matrix el¬ 
ement M{piP 2 9192) for Dtt Dtt scattering in 
Eqs. and i) can be written as 

T = Y\ f 

X { 2 tt)'^S'^{pi +P2-qi- 92) 'D'^M'D', ( 21 ) 


where 


M 

n 

!! 


{aDo{pi) aD+{pi)), = {a'^jyoiqi) a^+(gi)), 

ff^o(9i92;piP2)Tr(nn) i 


+Ma{qiq2-,PiP2)^^ + Me(gi92;piP2)nn, 
( a^+(P 2 ) \ 

y a-K-(j>2) -^a^o{p2) j ’ 


a 

71 

1 


-fe) -^alo{q 2 ) 


( 22 ) 


and ttM, Ojvf ^^e annihilation and creation operators for 
the meson M, respectively. Eq. (H is the most general 
isospin invariant operator for D{pi)T:(jp 2 ) D'{qi)Tr'{q 2 ) 
scattering. The T-operator is defined such that the famil¬ 
iar relation of T matrix (used in Eq. (H)) and amplitude 
M (used in Eq. (||)) can be reproduced: 


( 9 ig 2 |T|piP 2 ) = {‘ 2 TT)^S'^{qi + q2 - Pi - P2) 

xM{qiq2;piP2)- (23) 


, we have 


M{D+tt- D+TT-) 
M{D+n- D°n°) 


D°Tr°) 


Mq -\- Me, 

Mq + Ma , 

[Ma — Me] 

M(£)°7r° ^ D+tt-), 

Mo + ^[Ma+Me], (24) 


where Mi stand for Mi{qiq 2 -,P 1 P 2 )■ We can now make 
use of Eq. (||) to obtain ImM. For example. 


2 ImAjy0j^- = 21m M{B 
f (fqi (Pq2 


(27r)32Ti { 2 ttY 2 E 2 
xM*{D°Tr- D°tt-) Ado^- 
^ (Pqi (Pq 2 


- ^ 1)%-) 

{ 2 ti)'^ 5 '^{pi +P2- qi- q2) 


- J {2,m ( 2 ,)= 2 B. i (P. + B - 9 . - ®) 

X [Mo{qiq 2 ;piP 2 ) + M*{qiq 2 ;piP 2 )] Ajyo^- 
= {r*+rl)ADO^-, (25) 


where 


(Pqi (Pq 2 
(27r)32Ti { 2 ttY 2 E 2 


x{ 2 TTf 5 '^{pi +P2-qi- q2)M*{qiq2-,piP2), 


(26) 


for i = 0 , a, e and pi + P 2 = Pb- 

Since the Dtt system from B decay is S-wave, Aut^ is 
independent of the final state momentum. It can hence 
be factored out from the integration. Thus, Eq. ( [^ ) 
projects out the S-wave Dtt rescattering amplitude. Sim¬ 
ilar expressions can be obtained for ImAj^+^- Qo^o. By 
comparing with Eqs. (||) and (^, we find 


Eq. (II is unfamiliar as it is expressed directly in cre¬ 
ation and annihilation operators. It becomes more fa¬ 
miliar when expressed in fields, where we separate out 
creation and annihilation parts. For example, H and H 
correspond to the annihilation and creation parts of 



respectively. The SU(2) transformations of H and H can 
be recognized, since they do not mix creation and annihi¬ 
lation parts. One can find examples of using creation and 
annihilation operators on the studies of tttt, 7r-nucleon 
scattering in, for example. Ref. Q. 

It is important to note that the D{pi)tt{p 2 ) —*■ 
D'{qi)7T'{q 2 ) scattering amplitudes AI{Dtt D'tt') 
can be decomposed into the independent amplitudes 
Mo,aAqiq 2 ^PiP 2 )- For example, by using Eqs. (iD and 


ImAjjo^- ^ 

1 

0 


Im Ad+^- 

1 

+ 

II 

(27) 

ImAao^o , 

' V Ti-DOn° j 



rt=f '0 ' rS+r-: I, (28) 

0 -^(Ea-rl) + ^{r*a + r*e) 


with 5 = 1-1- iT. Eq. (^ ) is consistent with Eq. (|l^ ) 
through the identity 2Im5^/^ = for symmetric 

S. We also note that T can be diagonalized by using 
T = C/'^Tdiag U, where U is given in Eq. (|l^), giving 


7(iiag = diag(ro -I- Ve, Tq + re, cq -t- “ ’’e))- (29) 


The unitary scattering matrix 5 = 1 -f iT can 

be solved by identifying the elements in Tjiag with 
2 sin(angle) exp(i angle)’s, where we note that 1 -|- 
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FIG. 4: Pictorial representation of (a) (charge exchange), 
(b) Ta (annihilation) and (c) ro (singlet exchange). 


i 2 sin(angle) exp(f angle) = exp(f 2angle). One can now 
reproduce Eq. J^ ) by taking 

ro + re = 2 sin S 3/2 e®'^ 3/2 ^ 

ro + ^i^ra-re) = 2sin(5i/2 (30) 

We give a pictorial representation of Ve, ra, ro in Fig. 
they correspond to charge exchange, annihilation, and 
flavor singlet exchange rescatterings, respectively. Since 
the quark model is a representation of flavor SU(2) 
(SU(3)) group, it should be able to reproduce the struc¬ 
ture of T, which follows from symmetry argument. In 
particular the coefficients of in Eq. ( |^ ) can be repro¬ 
duced easily using this pictorial approach by matching 
the flavor wave function coefficients. For example, we 
have {va — re)l'/2 for a Z?+ 7 r“ to rescattering. We 

see from the first diagram of Fig. 4 that the exchange 
rescattering {je) projects out the dd component of 7 r° in 
the right hand side. From our convention in Eq. (|^), 
this give a — 1/-\/2 factor from the wave function of 7 r°. 
Similarly the second diagram of Fig. 4 projects out the 
uu component of 7 r° and give ra/ consequently. These 
diagrams also provide further information. For example 
in the second diagram as we go beyond SU(2), it is easy 
to see that the annihilation rescattering (r^) is responsi¬ 
ble for the D/'K~ rescattering, since there is 

no s quark before rescattering. 


IV. EXTENSION FROM SU(2) TO SU(3) 

A. DP Modes 

We now generalize the B Dtt case to SU(3) related 
modes in the final state, such as D°r] or DfK~. 

We stress that we apply SU(3) symmetry only towards 
final state rescattering rather than to the whole decay 
process. It is thus different from the usual application 
of SU(3) in B decays |Q, where one decomposes the 
B meson weak decay amplitudes, including the effective 
Hamiltonian itself, into different SU(3) pieces, and try to 
relate different modes (oftentimes including Bg decay). 
As argued earlier, SU(3) should be a good symmetry for 
energetic FSI rescattering, which is the case of interest. 
As we will see this approach give identical results with 
the SU(3) decomposition approach shown in Appendix 
A. 

It is st raig htforward to follow the steps through 
Eqs. (H) -(1^. Eq. (ill) remains unchanged, i.e. 

TT /■ d^p^ d^qj 

JX2' (2^)'2F;a2^)32F, 

X (27r)'‘<5^(pi +p2-qi- 92 ) D^MD', (31) 
but now the multiplets are extended 

pt = {aooipi) ao+ipi) a^+(pi)), 

= {a^jjoiqi) al,+ (gi) a|5+(9i)), 

M = Mo(qig 2 ;piP 2 ) Tr(nn) I 

+Ma{qiq2;PiP2)Bn + Me{qiq2;piP2)BIl, 


n(p 2 ) = 
/ 



/ 0 - 7^0 1 ^TJg 

\/2 Ve 

^-71“ 

Q'7r+ 

_|_ ^V8 

V6 

aK+ \ 

Qro 

(P 2 ), 



\ O.K- 

a-j^o 




n(g2) = 






( “lo , < 
\/2 Ve 

1 

+ 





®l+ 

4o 

( 92 ). 

(32) 


1 





Note that this operator can rescatter D^tt , 

into 


the desired states D'^rjs, Df K . 

The physical ry, rj' mesons are defined through 

f V \ _ f COS1? -sm§ \ f r]s\ 
yq' J sint? COST? / V / 

where the mixing angle 'd = —15.4° |^. In principle, 
we should also include D^r]i in the rescattering process. 
The additional terms can be obtained by replacing H in 
Eq. (^ ) by n -I-1 arix{q 2 )/Vi (and similarly for H), and 
labeling the qi related matrix elements by M[(ciiq 2 ] P 1 P 2 )■ 
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Knowing that the U^(l) symmetry is broken by anomaly 
and rji is not a Goldstone boson, M[ are not identical to 
Mi. The number of parameters would therefore double, 
but experimental measurements are still scarce. On the 
other hand, we note that the mixing angle i? is quite 
small, so we approximate rj hy rjs. Thus, we concentrate 
on the rescattering process involving octet pseudoscalar 
mesons only, as a step beyond the elastic FSI discussed 
in the previous section. In this way, as already shown, 
one again has just three independent amplitudes. 

Besides Eq. (|^) we now have 


M{D+Tr- D°r]8) 
M{D+Tr- D+K-) 

M(£)°7r° ^ 

DfK-) 
M{D %8 ^ D°r] 8 ) 
M{D%8 ^ DtK-) 
M{DtK- ^ DtK-) 


V6 




Ma: 

1 


{Ma + Me)^ 


2 V 3 

Mq + -^{Ma + Me), 
0 


{Ma - 2Me), 


1 

71 

Mq + Ma, 


(34) 


where Mi stands for Mi{qiq 2 ; P 1 P 2 ) ■ Amplitudes for other 
related modes can be obtained by noting that M{ab 
cd) = M{cd ab) in our case. 

With Eqs. (p4|), and similar extensions of Eqs. ( p^ ) 
and (p^, we extend Eqs. (27) and (p|) to 2ImA = T'A 
in the D^'k~ , D^irdD^rj^ and Df K~ basis with 


T = rol + 


/ 

0 

0 

0 

0 


'a+^e 

Ta 


0 

Ta—Te 

0 

T'a+T’e 

0 

\ 

d2 

TaA-Te. 


Ta 


2 

ra+re 

2 V 3 

ra+T’e 

V2 

ra—2re 


2\/3 

Ta 

6 

Va — 27*6 

V6 

) 

V2 

V6 

a 


(35) 


where T can also be easily obtained by the pictorial ap¬ 
proach, as we explain in the end of the previous section. 
Note that T can be diagonalized as 


"^diag = ro 1 -I- diag 





.(36) 


Following similar procedure in the previous section, 
5 = 1 -I- iT is obtained by the (physical) substitution 

(l + iro)e-2*^3/2 = + 

ira = i(-l-2e2*'5'-f3e^*®), (37) 

8 


where ro + re = 2 sin ^ 3 ji e®'^ 3/2 3,3 in Eq. ( po|) , but 
(5 = ( 5 i /2 — (^ 3/2 is now extended to two physical phase 
differences 8 ' and 9. Note that D, Dg and tt, K, rj trans¬ 
form as 3 and 8 under SU(3). As shown in Appendix A, 
we can identify above phases as 

^ 3/2 = S = 60 — djg, 9 = 6 -^— (38) 


By analogy to the previous section, the solution of 
2 Im A = A is 


/ \ 

■^D+-n- 

V ^DtK- / 


= 51/2 




A 

A 

A 

[< 


jP+TT- 

D°ris 


D+K- 


\ 

/ 


(39) 


where A-^s are factorization amplitudes. The matrix Sif 2 
can be obtained by reducing phases in S by half. For 
later purpose, we give the explicit expression of (or 
equivalently S with trivial modification of phases). 


^5l/2g-i<53/2 

/ 1 0 0 0 0 \ 


0 


g^(-5 + 2e*^'+3e*«) 


-i(l-h2e*'^' -3e*®) 


0 

^(-5 + 2e*^'+3e*®) 

jh(ll + 2e*'^' -f 3e*®) 



.(40) 

0 



-^i9 + 6 d^' +e^^) 

i V^(-3-f 2e*'5' -h e*®) 


VO 

-|(l-k2e*^' -3e*®) 



i(3-k2e*^'-k3e*®) yi 



Just as in Eqs. (^ and (p^), 5^/^ of Eq. ( p0| ) has an 
overall phase. Only phase differences affect decay rates. 
An overall sign change of the phases also leaves rates 
unchanged. 


Note that charge conservation and unitarity imply, 
|Ado^- P = \A^jjo.„- 7 

\^D+'k- P + P + |Ado^P -I- |A^+^_ p 

= + |A{jo^oP + (41) 
























Since the amplitudes for color suppressed modes are 
small at the factorization level, they will be fed by the 
color allowed amplitude As a consequence, the 

rate will be reduced from its factorization result. 
Compared with the elastic case, we now have additional, 
but slight, leakage of D^tt~ into D^rj and Df K~ modes. 
Because the measured rates of color suppressed modes 
are still small as compared to the color allowed modes, 
as far as the Dir system is concerned, the results do not 
deviate too much from the previous section. 

On the other hand, since the factorization amplitudes 
satisfy the isospin triangular relation, Eq. (p^), one 
can show that, with ESI and for any value of and 

the rescattered amplitudes also satisfy the 
relation. As noted earlier, the isospin relation should 
hold whether ESI is active or not. To demonstrate this, 
it is instructive to express the ESI of Eq. (|4^ ) in the 
isospin basis. In the D^tt~, D^rjs and 

Df K~ basis, the T matrix, and similarly take a 

block diagonal form 


i^lock — 


/ 

I’e 

0 

0 

0 

0 

\ 


0 

re 

0 

0 

0 



0 

0 

3ra—re 

2 

ra+re 

2 

V 2 ^ 



0 

n 

Ta+re 

ra+r-e 

ra—2re 




2 

6 



V 

0 

0 

\l 2 ^ 

ra—2re 

V6 

ra 

/ 


(42) 


elastic” rescatterings 


g-*53/2 150,,^ = 


p-iS3/2 

^ l‘^blockl(D7i-)i/2,r>4if- 

g-*<53/2 

p-iS3/2 


('^block) D0r,s,DtK- 


a-i^3/2 'l , , 

l^blocklD+if- n+/f- 


l(l + 6e*49e*®), 

lo 

l(9 + 6e*4e*®), 

16 

^ (-3 + 2e*'*'+e*®), 


8v/2 


= ^(3 + 2e*'^+3e*®). (45) 


between the / = Iz = 1/2 decay final states 

(Il7r)i/2) and D^K~, which is a reasonable ex¬ 

tension beyond the elastic rescattering discussed in 
the previous section. The elastic case corresponds to 

e~'‘^^^HK(ock)(D7r),/2.(D7r),/2 = and setting the rest 
of Eq. (^ ) to zero. Since the factorization amplitudes 
for D°77 and K~ are small, their ESI contribution to 
the D-k system will be suppressed. Therefore, 


2^1/2 


^1/2 



(1 -k 6e*^+ 9e*®) 


arg (1 -k 6e*'*'+9e*®) 


(46) 


are good estimates. The geometric meaning of Eq. (|4q ) 
is given in Appendix B. 


by 7),iock = OTO’’’, where 


B. D*P and DV Modes 


t/ 0 0 

0=1 0 10 
0 0 1 


(43) 


and U is given in Eq. (18). Truncating to the first 3x3 
sub-matrix, T is diagonal and one reproduces the SU(2) 
case of Eq. 


1 /2 

In this basis, with block diagonalized = 

Eq. (p^ becomes 


The formalism can be applied to final states involving 
pseudoscalar and vector mesons (PV), with only slight 
modifications. 

For D*{pi;X)P{p 2 ) D*{qi\\')P'{q 2 ) rescatter¬ 

ing, where is the polarization index and P^''> de¬ 
notes a pseudoscalar meson, we replace anipi) and 
Mi{qiq 2 ;PiP 2 ) by aD.(pi,A) and Mi{qiq 2 , X';piP 2 , X) 
etc., respectively, where 

(9192, A'|T|piP2, A) 

= {2Tr)'^S'^{pi -kp2 - 9i - q2)M{qiq2,X'-,piP2,X). (47) 


A _ cl/2 ./ 
^ — ‘^block ^ > 


(44) 


The B D*P amplitude is expressed as 
M{B ^ D*P)=el-pBAD^P, 


(48) 


where the first two diagonal elements of 


g*<53/2^ while the remaining lower block governs the “in- 


where e\ is the polarization vector and A]j*p is a Lorentz 
scalar that is independent of e and the angle between the 
3-momenta of D* and P. 
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Take B D*^tt for example. Eq. 


becomes 


= E 


2Im[(e^ ■Pb)Ad.o^-] 

'■ (Pqi (Pq 2 


{2T:f2Ei (27r)32E2 


{2t:Y5^{pi +P2-qi-q2) 


X [M^{qiq2,\']piP2,X) + Ml{qiq2,\']piP2 A)] 

xsy ■ pb (49) 


By choosing a real basis for e and noting that ■ 

e\)^ = TOg/m^.o, where Pcm is the momentum of 
in the center of mass frame, we obtain 

21mAB.o^- ={r;+r;)AB.o^-, (50) 

with 

* = V f 

~ {2TTf2Ei [2 ttY2E2 

y.i2.f5Hp^+P2-q^-q2)^^ 
xex ■ pb M*{qiq2X';piP2, X) £y ’Pb, (51) 

which projects out the P-wave D*P scattering amplitude. 
Equations after Eq. (|^) from previous section can be 
carried over by replacing Ajjp Ap^p. 

The generalization to DV decay modes is again 
straightforward, except that wi-wg mixing cannot be ne¬ 
glected. The physical mesons should be 


CO = 



(52) 


Unlike the pseudoscalar case, where ignoring rji couplings 
can be partially justified, coi in principle should be on 
similar footing as wg- 

We replace tt, K and rjs in Eq. ( ^^ by p, K* and wg. 
By including coi we will have two more terms in A4 of 
the corresponding r matrix namely, al^^{q 2 )ac^^{p 2 )[Mo + 
(Ma + Me)/3] and (pz) B + B at,^ (g 2 )](Ma + Me)/V3. 

These two terms correspond to coi to coi and coi to octet 
rescattering. These Mi, Mi will reduce to Mi under 
nonet symmetry. 

The T matrix can be obtained as before. In isospin 
basis with the (T>p) 3 / 2 , (^p)i/ 2 , D^co^, DfK*- 

and D^coi, we have 


^lock = diag(ro, tq, tq, Tq, ro, tq) 


+ 



0 

0 

Tg+re 

2 

Tg+re 

6 

ra—2re 

V6 

r-g-l-r'e 

3\/2 


2 'a 

Tg —2re 

V6 

Ta 

rg+fc 

V3 


0 

0 

rg+fg 

V 2 

rg+Te 

3 V 2 

fg+fe 

rg+Pe 


.(53) 


It can also be easily obtained by the pictorial approach 
in the coi, wg basis. Tbiock is identical to Eq. (|4^ ) except 
the additional sixth row and column. The solution is, 
ro,e as shown in Eq. (^^ and 




1 


= ^( 3 ^^ 336 “^**"/" - 2e2*'5' - 1 ), 


ii.ra + Ve) = 


, ra+re^ 

^[ro +- 5 -) 


^3'3' ~ 


(54) 


where U is a, two by two symmetric unitary matrix that 
mixes 3 and 3 by rescattering as shown in Appendix 
A. Note that we need two phases and one mixing angle 
to specify U, resulting four parameters (with an overall 
phase factored out) to describe the rescattering matrix. 
There are too many parameters and experiment mea¬ 
surements are still scare. However, there is no U^(l) 
anomaly in the vector sector, and we expect U(3) rather 
than SU(3) symmetry. Therefore we consider rescatter¬ 
ing among and the vector nonet as a first step 

beyond the simple elastic ESI case. In this case we iden¬ 
tify fi and Ti as r^. 

It turns out that U(3) symmetry allows only either 
charge exchange or annihilation ESI, but not both. The 
two distinct solutions for and U are 


Solution 1: 


(1 -k iro) 

ive 

iva 

Z^33,e-2*^3/2 

Z^3,3,e"2*^3/2 

Solution 2 : 

(l-kzr^'^)e“(^^*^3/2 

g—(2)*53/2 
Ij-(') g-(2)ii53/2 


i(l + e^*^'), 

^(1-e^*''), 

0 , 

1(1+ 26^*''), 

^(1-e^*''), 

1(2 Te^^^'), (55) 


1 , 

0 . 


l(-l + e(2)-), 


3 
1 

9 

2 V 2 


(l + Se^*'^), 

(-l + e^-). 


-(S + e^-), 


(56) 


where the former (latter) does not have annihilation (ex¬ 
change) contribution. Notice that while 6 ' = Sq — (Jyg is 
analogous to the D^Ap counterpart, a is from hi and is 
not equivalent to 6 . 
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To understand why raVe = 0, we shown the T matrix 
in the D^p~, p~, D^p^, D^uj, Df K*~ and basis. 


Tv = tqI 


( 

0 
0 
0 
0 
V 0 


Ta 

Vg—re 

V2 

r-g+re 

72 


0 

Vg —Te 

V2 

Tg+re 

2 

Tg+re 


72 

0 


0 

rg+re 

2 

r-g+re 


72 

0 


0 

ra 

Vg 

72 

Tg 

72 

Ta 


0 \ 
0 
0 
0 
Te 

0 y 


(57) 


Note that Ty' = 0 for i = 2,3, 4, so can only rescat¬ 
ter with itself and D^K*~. This can be easily checked 
by using the pictorial method as shown in Fig. 4. 

We define Ty via = 1-1- iTy. It is easy to show 

that Tv — 2Ty iTy. should also satisfy the U(3) 
symmetry as well since it is generated by the same dy¬ 
namics (or Hamiltonian) as Sv- Since the construction 
of Tv is based on symmetry, we expect Ty to have the 
same structure, or simply with ro,a,e replaced by Cq g. 
It is then easy to show that 




(58) 


for j = 2,3,4. Since these elements of Tv are zero, we 
must have r' = 0 which implies ra rg = 0. 

With these r', which are in Eqs. (|55|), (js^ ) with 
phases reduced by half, we obtain the decay amplitudes 
by applying A = Af. 


V. NUMERICAL ANALYSIS 

A. R —> DP and D*P Modes 

For B Dtt modes, we start from Eq. ([l^), where 
the factorization amplitudes are decomposed into color- 
allowed external W-emission (T), color-suppressed inter¬ 
nal IT-emission (C), and W-exchange amplitude (E). 
They are given by 

T = ^ycbV*dai {m% - 

C = (59) 

Since the annihilation form factor Fq 

pected to be suppressed at and 02 is small, 

the amplitude E is neglected. For we have 

= ^ycbV*aa2{ml - to2)/dFo^’’(to^), (60) 

where pi-ps niixing effect is included via FQ^{rrijj) = 
cos d Fq'^^ — sin -d F^^^. We use experimentally measured 


TABLE I: Form factors in LF and NS form-factor models 
where Af‘^{q^) = AfP{q^)lV2 and V^‘^{q^) = V^'’{q^)ly/2. 



LF (NS) 


LF (NS) 

Fo^^iml) 

0.29 

(0.27) 

A^^\mi) 

0.73 

(0.64) 

Pf-(ml,.) 

0.34 

(0.32) 

A^»{ml) 

0.35 

(0.31) 


0.16 

(0.15) 

Af^iml,) 

0.23 

(0.28) 


0.19 

(0.18) 

A^'^iml,) 

0.22 

(0.31) 

r-iBD \ 

0.70 

(0.63) 


0.30 

(0.32) 

-^0 \^k) 

0.70 

(0.64) 

A BD* /^2 \ 

Aq \^k) 

0.74 

(0.65) 

r-iBD \ 

-^0 (^K*) 

0.71 

(0.65) 

aBK* (2 \ 

^0 \^d) 

0.40 

(0.35) 

TtBK \ 

^0 \^d) 

0.42 

(0.31) 

T?BK(^2 \ 

Fi 

0.44 

(0.36) 


masses in ^do^o and These am¬ 

plitudes are real in our phase convention. 

For B D*Tr, D*^p, we have M{B D*P) = 
{e ■ pb) Ad*p [Eq. (^)]. Analogous to Eq. ( 0 ) , the fac¬ 
torization amplitudes are decomposed into 

T = ^VcbV:^aif.21710^ A^^irnl), 

C = ^W6lC^a2/zt*2mz3.Ff"(ml,), (61) 


where E is again neglected. Eor B^ D*^p we have 

4 . 0 , = ^V,,V:^a2fD^2mD-F^\ml). (62) 

Starting from ESI redistributes these sources 

into the amplitudes Ajj{„)p, and we obtain the corre¬ 
sponding rates by, 

r(;B^PP) = |App|Vm/(8^m|), 

r(:B^UP) = \Avp\^plJ{ST:ml), (63) 


where 


2 mppcm = Y ["i| - (wi + m 2 Y][m?B ~ “ "12)2]. 

(64) 

Note that we factor out e ■ pp in the definition of the 
amplitude Avp, so our expression for F(P ^ VP) is 
slightly different from that in Ref. [ p^ . 

In our numerical study, we fix Vud = 0.9749, Vus = 
0.2225, Vcb = 0.04, and use the decay constants /tt = 133 
MeV, fiiV) = 158 (214) MeV, /p(.) = 200 (230) MeV and 
fp = 210 MeV. Masses and lifetimes are taken from Par¬ 
ticle Data Group (PDG) ||^. We consider two form factor 
models: the relativistic light-front (LF) quark model jl^ 
and Neubert-Stech (NS) model |^. The relevant values 
are listed in Table |. We use the color suppressed branch¬ 
ing ratios of Belle Q, except for modes where 

we combine with the latest CLEG numbers ||] . Eor other 
modes we use PDG values [^. Since the charged 
mode does not rescatter to other modes, we normalize 
all modes to = [5.3 (4.6) ± 0.5 (0.4)] x lO’^ 
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TABLE II: The best fits in the SU(3) FSI picture. The 
subscript indicates DP or D*P modes. Form factor model 
dependence is less than a couple of %. 



Fitlpp 

Fit2DP 

Fito'P 


0.20 

0.27 

0.21 

s' 

47.8° 

17.1° 

55.7° 

0 

24.8° 

-52.7° 

18.2° 

(1 -l-iro)e"^“''3/2 

0.45 -K 0.50i 

0.91 -h0.28i 

0.32 -h 0.47i 


0.55 - 0.50i 

0.09 - 0.28i 

0.68 - 0.47i 


0.14-H0.04i 

-0.43 - 0.50i 

0.27- O.Oli 

IA 1 / 2 /(■\/2A3/2)| 

0.75 

0.65 

0.71 

\Sl/2 — ^ 3 / 2 ! 

30.2° 

26.2° 

28.3° 


TABLE III: The branching ratios of various D^*^P modes in 
10““^ units. The second and third columns compare experi¬ 
ment with factorization model. The last two columns give the 
best fit results with FSI parameters of Table ^ 


Mode S(xl0-^) fac^^(N®) Fitlop Fitp.p 

- - 0.20 0.21 


D+tv- 

30 ±4 

35.7 

(35.5) 

32.2 

- 

D°7V° 

2.9 ±0.5 

0.57 

(0.58) 

2.93 

- 

D°ri 

1.4 /gi±0.3 

0.33 

(0.34) 

1.43 

- 

DtK- 

< 2.4 (0.7) 


0 

0.03 

- 

D*+-k- 

27.6 ±2.1 

29.8 

(29.0) 

- 

26.3 

D*°tv° 

2.5 ±0.7 

0.60 

(0.69) 

- 

2.44 

D*°ri 

2.0 ±0-4 

0.34 

(0.39) 

- 

1.83 

D*+K- 

< 1.7 


0 

- 

0.16 


1^. We then perform a fit to the ratios of branching 
ratios B{D<^*'>°Tr°)/BiD<^*^°TT-) 

and B{D^*^°r])/B{D<^*>°TT-). 


The use of ratios reduces model dependence on form 
factors, and is sensitive only to ai/a\. Our numerical 
results for rescattering phases in LF and NS form factor 
models never differ by more than a few %. With |a 2 | = 
0.26±0.02 from fit to B ^ J/ipK data and the range 
of oi ~ 1-1.1, 02 ^ 0.2-0.3 from various modes |^, we 
shall adopt 02/01 = 0.25 in subsequent discussion. We 
find that a larger 02/01 is preferred for the DP modes, 
but the converse is true for D*P modes. However, for 
02/01 ~ 0.25, the x^s of best fits to DP and D*P modes 
are both quite small. 

The best fits for FSI phase differences 5' and 9 (or 
alternatively the rescattering parameters r^) are given in 
Table ||. We do not list form factor model dependence 
since it shows up often only at third decimal place. We 
find two fits for the DP case, but only one fit for the D*P 
modes. For DP modes, the set that we call “Fitl” is 
similar to the D*P modes, i.e. 6 ' ~ ±50°, 6 ~ ±20°. We 
find that the quark exchange strength |re| is larger than 
the annihilation strength |ra| in this case. As illustrated 
in Appendix B, the large phase ~ 50° arises because of 
sizable strength of both and While S' and 

6 effects are of similar sign for the former, for the latter 
they counteract, and a large <5' phase is needed. 

For “Fit2” of the DP case, we have S' ~ ±20°, 9 ~ 
=F50°, implying that \ra\ > \re\- As we will see later, 
this fit is ruled out by the Df K~ bound. We 

mention a curious point about our “Pomeron” related 
effect, i.e. |l±iro| = 0.67, 0.95 and 0.57, respectively, for 
Fitl/jp, Fit2pip and Fitp.p. The first and the last are 
remarkably consistent with the estimates of >- 7 rir ~ 
0.58 @, 0.68 and ~ 0.76 § at Vs = mp. In 

contrast, Fit2pip, which is already ruled out by data, is 
not quite consistent. 


In Table || we show |Ai/2/(V2A3/2)| and \Si /2 - < 53 / 2 ! 
obtained by usiM the fitted strong phases. The compar¬ 
ison with Eq. dig ) will be discussed in the next section. 

We summarize the predicted rates of various modes 
in Table IB. The branching ratios are obtained by 
multiplying the fitted ratios of branching ratios by the 


measured central value of B{D^*'>^Tr~). The results for 
T)(*)+7r“, fit the data well, as expected 

from the small x^. It is worthy to note that rates given 
in Table Bl satisfy the unitarity condition of Eq. (jj), 
i.e., the sum of branching ratios before and after FSI is 
equal as shown in Eq. (|l]), up to small phase space cor¬ 
rections. This clearly shows that the color suppressed 
modes are fed from the D^*'>~^tt~ mode. 

At the amplitude level, one reads from Eq. ( ^ ) that 
Il(*)°7r°, D^*^^ri and Di*''^K~ receive the (r^ — re)l'/2, 
(^a ± ?"e)/'\/6 and rescatterings, respectively, from 
£)(*)+ 7 r“. Indeed, by using phases shown in Table || 
in Eq. (|40|), the FSI contribution can be estimated by 
using |Ap(.)o,ro/AD(.)+^-| ~ and sim¬ 

ilarly for |Ap)(.)o,^| and |A^(.)+^_|. The FSI contribu¬ 
tion from D^*'^^Tr~ alone provides 70%-80% of measured 
7r°, D^*^^ri rates, with the remainder coming from 
02 and interference terms. For the Ds K mode, the 
FSI contribution from is small due to the small¬ 

ness of Ta but it is still three times larger than that shown 
in Table IB, Because of this smallness, the FSI rescat¬ 


tering from the color suppressed modes (due to the non¬ 
vanishing 02 ) through Te cannot be neglected, which re¬ 
duces the rate to that shown in the Table. 


We note that the factorized 02 contribution to color 
suppressed modes show some form factor model depen¬ 
dence, especially for D*P modes, but such model depen¬ 
dence for fit results are rather slight. The reduced form 
factor dependence is quite consistent with FSI rescat¬ 
tering dominance over factorized 02 amplitude, since 
d(*)+ TT is the common source. 


In “Fit2” wetindB{D+K-) = 4.65xIO"'* > B{D°n°), 
B{D^r]) due to \ra\ > \re\- This value is above the PDG 
bound of B{DfK~) < 2.4 x 10““^, and way above the 
recent Belle bound of 0.7 x 10“"^ [^, hence is ruled out. 
The results for Fit2 are therefore not shown in Table IB. 
On the other hand, since |ra| in the D*P case is not too 
small, the result of B{Dl~^K~) ~ 1.6 x 10“® may still be 
of interest. 
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B. FSI in DV modes 

For B Dp, Duj, we have M{B ^ DV) = e-ps Adv 
and can be decomposed into 


C = ^VcbV*aa2 fD^mpAgPiml,), (65) 


while E is again negligible. For 5° —> D'^uj one has 


02 fD2m^A^-{ml). ( 66 ) 

We set = 0 for the factorization ampli¬ 

tudes of B^ —> DfK*~, D'^cj) modes. We again normal¬ 
ize to D^p~ mode since its rate is unaffected by FSL We 
take the D^uj measurement from Ref. 0 , while the mea- 

-\)p- 

B{D^p^) are taken from PDG 

Because of the reduction to one phase difference for 
both Solution 1 and Solution 2 of Eqs. ( |^ ) and (p^), we 
are able to fit with just D'^p~ and D'^uj data. We find 6' 
(a) to be 18° (35°) for Solution 1 (Solution 2), as given 
in Table Since the FSI contributions to D^p^ and 

-Af 

^DOuj^ 


surements of D^ p , D^p modes and the upper limit of 


D^lo are mainly fed from H+p 
Eq. (p^) leads to 


and 24{,0p0 


^D°p°(ui) 


^KTr',) f 

/K ^D+n- + G 


V2 


DOpO(ij) 


^ ^D+p-T (1 + (67) 


where r'^, r'^ and are defined in Sec. IV.B. For So¬ 
lution 1 (Solution 2) the dominant FSI contributions to 
D^p^ and D^uj are the same in magnitude but oppo¬ 
site (same) in sign due to ri''^ 7^ 0. This implies 

different interference patterns. For Solution 1, we have 
B{D°p°) ~ B{D^uj). For Solution 2, since and 

^z)0aj same sign, 

(I -t = 1, 

Re ^ [cos cr — I] < 0, (68) 

O 


TABLE IV: The best fit phase difference for DV modes. 


Mode 

Solution 1 

Solution 2 

x"" 

1.20 

0.64 

5' 

18.1° 

- 

a 

- 

34.8° 

(1 -bfro)e"^“'3/2 

0.90 + 0.30i 

1 


0.10-0.30i 

0 

irae~'^^A/2 

0 

-0.22-b0.3H 


TABLE V: The branching ratios of various DV modes. The 
second and third columns compare experiment with factor¬ 
ization model. The last two columns give the best fit results 
with FSI parameters of Table s 


Mode 

B (xlO-^) 

facLF (NS) 

Solution 1 

Solution 2 

D+p- 

79± 18 

100.7 (101.2) 

98.2 

92.7 

D°uj 

1.8 ± 0.5 j 

0.67 (0.64) 

1.86 

1.92 

D°p° 

< 3.9 

0.67 (0.64) 

1.90 

3.37 

DfK*- 

< 9.9 

0 

0 

2.73 

D°(f) 

- 

0 

0 

0 


the FSI contribution always interferes destructively (con- 
tructively) with A^j^g^ {A^j^o^q). While B{D°p^) becomes 
larger, one would need large annihilation contribution to 
account for the observed D^uj data, which in turn gives 
rise to B{Df K*~) as large as 2.7 x 10“^. These patterns 
can be tested in the near future. 

On the other hand, 1?°^ only rescatters with £)+ K*~, 
as can be seen from Eqs. (p^). It does not pick up any 
FSI contribution since A^ . _ = 0, even if A^+j^,_ is 

nonvanishing as in the Solution 2 case. Observation of 
D^ip mode would imply some mechanism at the “source” 
level. Our fitted branching ratios and predictions for var¬ 
ious PV modes are given in Table 0 

We note that [711/2/(^2^3/2)1=0.84, \5i/2 - < 53 / 2 ! = 
12.7° from Solution I and |7li/2/(-\/2A3/2)|=0.8I, \5i/2 — 
< 53 / 2 ! = 18.5° from Solution 2. The phase angles are 
somewhat smaller than those for D^*'>'k modes. 


C. Predictions for D^*^K and DK* Modes 


Our FSI formulas can be applied readily to rescattering 
between the Cabibbo suppressed modes, D^*'>^K~ and 
D^A^K^, since they are contained in the formalism for 
DV)p. Following similar procedure as before, we have 


Vlm^D(*)o7fo / V / 


(69) 


where 


r = 


ro Te \ 
re ro ) 


(70) 


and “annihilation” is clearly impossible as can be seen 
from the pictorial approach. By using Eq. (^), we note 
that 5 = 1 -b iT is automatically unitary. Therefore, we 
obtain 


( Ajj(,)+ii- 

V Ajj{*)oxo 


gl/2 f ^dM + k- 
\ ^DMoko 


(71) 


where. 


= 


gi'53/2 


1 -b 1 - \ 

1 - 1 -b ) ’ 


2 


(72) 
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TABLE VI: The predicted branching ratios of D^*^K and 
DK* modes in 10“^ units. The second and third columns 


compare experiment with factorization model. The last three 
colnmns give the best fit results with ESI parameters of Table 


and Table IV 


Mode 

B (xlO-4) 


Fitlop 

Fit D*p 

Solution 1 

D'^K- 

4.2 ±0.7 

4.12 

- 

- 

- 

D+K- 

2.0 ±0.6 

2.61 

2.20 

- 

- 

D°K° 

- 

0.12 

0.53 

- 

- 

D*'^K- 

3.6 ± 1.0 

3.40 

- 

- 

- 

D*+K- 

2.0 ±0.5 

2.15 

- 

1.71 

- 

D*°K° 

- 

0.10 

- 

0.55 

- 

D'^K*- 

- 

7.0 

- 

- 

- 

D+K*- 

- 

5.11 

- 

- 

4.99 

D°K*° 

- 

0.09 

- 

- 

0.21 


which is consistance with those obtained in Appendix A, 
and the factorization amplitudes are 


^D+K- 


DOKO 

Af 

^D* + K- 


A^ _ 

D*0K0 


Gf 

V 2 

Gf 

Gf 

G 


VcbV*, ai {ml - 

VcbV*^ 02 {ml - TOxo)/r>^o^^(w|,o), 
VcbV*^ oi 2mD-fKAQ^{ml), 


= -^V,bV:,a22mD>fD*F^^{ml.o), (73) 


72 


One again has 
^ DMoko- 


where form factors are found in Table 
a triangle relation 

The 6 ' phase has already been fitted, and the predicted 
branching ratios for DK and D*K modes are given in 
Table VI. The second column is obtained by multiplying 
Belle measurements of B{D^*''^''^K~)/B{D^^’^f~) 
by PDG values of B{D‘^*'^+'°f-) [|._ 

It is interesting to note that IB*^K modes do not re¬ 
ceive the annihilation type FSI, and the wave 
function does not have the l /-\/2 factor as compared 
to 7 r°. For “FitlDp” and “Fito*?, Va is subdominant 
while I Tel is close to each other, hence we find that 
B{D°K°)/^D+K-) « 2B{D^Tr°)/B{D+F-) « 2 x 1/10 
and B{D*°K°) « B{D°K°). As noted, “Fit2Dp” is ruled 
out already by Df K~ bound. 

In a very similar fashion, we predict the rescattering 
of D'^K*~ into D^K*^ final state, which is given again 
in Table 0 for Solution 1 of DV case, where one again 
hasB{D°K*°)/B{D+K*-) « 2B{D°p°)/B{D+p-) « 2x 
1/50. For Solution 2, Cg = 0 and the result is the same 
as the second column for factorization. 

The DK* modes have yet to be observed. The fac¬ 
torization predictions for D^K*~ and D'^K*~ are about 
twice as large as D*'^K~ and D*~^K~, but the predicted 
branching ratio for B D^K*^ is less than half of 
B D*^K^ in Solution 1. This is because |re| (or the 
rescattering phase S') for DV modes are much smaller 
than for D*P modes. We note that in Solution 2 one 


would predict B D^K*'^ to occur at half the rate of 
the Solution 1 case, i.e. just the factorization 02 predic¬ 
tion. 

Our formalism therefore predicts a relatively sizable 
B{D<^*'>°K^) at - 0.5 X 10-4, and expect B{D°K*°) - 
0.2 X 10-4 (0.1 X 10-4 if Solution 2 is confirmed). We 
encourage Belle and BaBar to search for these modes. 


VI. DISCUSSION AND CONCLUSION 


We have stressed that the isospin relation of Eq. dil), 
which follows from Eq. (|^ (given more explicitly in 
Eqs. (P)-([ll|)), holds whether one has FSI or not. This 
was used to extract ^ 1 / 2/^372 directly from D^*^f data, 
as given in Eq. ([T^, which we reproduce here: 


1 

7 I 


^ 1/2 

^ 3/2 


TT only 


l<5| 


only 


0.71 ±0.11 (0.75 ±0.08), 

29° ± 6 ° (30° ± 7°), (74) 


where it is made clear that they are extracted from D^*'>ti 
data alone. On the other hand, we have given in Ta¬ 
ble H the values for |>li/ 2 /( 72 A 3 / 2 )| and \Si /2 - < 53 / 2 ! 
as obtained by using the fitted strong phases that takes 
into account D^ri^*F data, i.e. by using the Ajj(,)o^-, 
Aoi *)+^- and Aj^(*)o^o amplitudes of Eq. (Il9|) . They 
turn out to be not so different from Eq. ([74|). Let us 
understand why. 

We note that |xli/2/('\/2A3/2)| = 1 ± ©(A/mg) —> 1 
in the heavy b (and c) quark limit &0 , although me may 
not be heavy enough. The strong phase |(5i/2 —< 53 / 2 ! —^ 0 
in the absence of short and long distance rescattering. It 
is intructive to consider first the elastic D^*'>Tr DG) 

TT 

rescattering case. Noting that elastic rescattering does 
not change |A/| from its factorization value, for D^'lp 
modes we have, roughly speaking, 


I-41/2I l 4 / 2 l | 2 r-C| | 2 -gf| 

72IA3/2I 72IA3//2I 2 |T±Cr2|l±f^|’ 


which deviates from 1 due to a nonzero 02 , which is a 
nonfactorizable effect. It was a happy coincidence, before 
the measurement of color-suppressed modes, that taking 
02/01 = 0.25 and real could account for [Q both B~ 
D^Tr~ and B J/ipK^) rates. It should be stressed that 
the sizable value of |a 2 /ai| can be viewed as determined 
this way from data that gives |Ai/ 2 /(-\/ 2 A 3 / 2 )| — 0.7. 

The impact of the new experimental measurement of 
D°tt° is that one is now able to determine the strong 
phase difference, \ 5 i /2 ~ ^ 3 / 2 \ 30°j which is not quite 
small. With this one has two ways to proceed. 

As mentioned in the Introduction, Refs. [D, 0 con¬ 
tinue to employ factorization formulas to Il^*™ 7 r"hence 
make |a 2 | larger by roughly a factor of two. To main¬ 
tain dG)^ TT and generate 61/2 ~ ^ 3/21 ons capitalizes on 
Eq. ([^ by dropping the absolute value condition. That 
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is, one resorts to a complex a^ja\ itself. In this way one 
finds \a 2 \ 0.4-0.5 and arga 2 ^ 60° could account for 

B Dtt data, at the cost that |a 2 | is twice as large 
as from modes. Ref. |0] argues further that, 

while factorization no longer holds, the trend of larger 
and complex 02 is expected from QCD factorization g. 

Our critique is the following. First, this “ 02 ” ap¬ 
proach is process dependent, and predictiveness is lost. 
Although |a 2 | ^ 0.4-0.5 could account for the strength of 
observed D^*^^r] and uj modes, it seems coincidental, 
with 1 02 1 varying by ~ 20-30% among these and D^*'>^Tr^ 
modes, while we know that |a 2 | ~ 0 . 2 - 0 .3 for 
modes. Second, it is the need to maintain rate 

that a sizable phase to 02 is invoked, although previ¬ 
ously a smaller and real 02 gave a pleasant, consistent 
picture. However, from Eq. ( [7^ we know that 02 reduces 
|Ai/ 2 /('\/ 2 A 3 / 2 )| hence it represents inelastic effects 
This is the reason why it is not quite calculable. In this 
sense, however, arg 02 ^ 60° is not reasonable since one 
expects strong cancellations among numerous inelastic 
channels A statistical model suggests the typical 

phase to be ~ 20° j|]. Third, we stress that generat¬ 
ing Eq. (^) by the phase and strength of 02 holds only 
when one drops the absolute value requirement from Eq. 
(^), i.e. ignoring elastic ESI, as is common practice in 
QCD factorization |^. Such ESI effects are 0{as{mi,)) 
suppressed, or IfiriQ suppressed. For the former, clearly 
as{mb) ~ 0.2 ~ 10° in radians. For the latter, one has 
the real problem that rric may not be heavy enough. 

The approach advocated in this paper is via quasi¬ 
elastic FSI rescattering. Let’s make a point by point 
comparison with the “process dependent 02 ” approach. 
First, process independence is not so easily conceded. In 
particular, we maintain 02/01 = 0.25 and real. Thus, 
the proximity of Eq. with Eq. reflects a phe¬ 
nomenological tuning done several years ago. The phi¬ 
losophy is that, while we agree that 02 should be process 
dependent and in principle complex, we take the above 
“tuned” value as no mere accident. That is. Nature could 
have revealed to us long ago that 02 is strongly process 
dependent. Since keeping both 02 as parameter and al¬ 
lowing for FSI phases cannot lead one afar, we opt to 
keep 02 !ai fixed as done in 0. The new experimental 
measurement of and modes, rather than 

giving “process dependence” irritation, can be incorpo¬ 
rated nicely by enlarging the scope of FSI from elastic 
SU(2) to quasi-elastic SU(3) symmetry. This stretching 
of “elasticity”, together with maintaining process inde¬ 
pendence to good degree, make our approach suitably 
redictive. Second, by leaving 02/01 as done before in 
, one enjoys the success with D^tt~ and J/ipK. 
The strength of 02 is smaller, hence the acuteness of in¬ 
elasticity is not as severe as the “ 02 ” approach. A strong 
phase of order 50° does emerge, but this is interpretted 
as due to having D^* )° 7 r° and rj both sizable (see 
Appendix B), and has completely different origins from 
the need for arg 02 ~ 60° in the large “ 02 ” approach. 
Finally, in comparison with the strong assumption of re¬ 


moving the absolute value condition from Eq. (|75|), we 
took advantage of FSI approach to expand a previously 
commonly known “folklore” on elastic FSI rescattering. 

We can now comment on comparison of Eq. (|7^ ) 
with values in Table II. With 02/01 = 0.25, we find 
|Ai/ 2 /( 72 A 3 / 2 )| = \A{^JiV24^^)\ = 0.77 (0.75) for 
LF and 0.77 (0.73) for NS form factors. Assuming elas¬ 
tic rescattering and using the formulas of Sec. Ill, the 
strength of amplitudes cannot change, and we find \5i/2 — 
53 / 2 ! = 30.6° (29.9°) for LF and 30.6°(29.2°) for NS form 
factors. Note that the form factor dependence is very 
weak. With rescattering among SU(3) multiplets, 
can now feed other color suppressed modes via rescatter¬ 
ing (Eq. (H)). The TT mode still can not rescatter 
to other modes, so IA 3 / 2 I = ^ 3^2 = \Ado-k-\/ Thus, 

from unitarity we expect IA 1 / 2 /A 3 / 2 I < 1 ^ 1 / 2 /^ 3721 
the presence of quasi-elastic FSI. We see from Table || 
that IA 1 / 2 /A 3 / 2 I is reduced from the factorization results 
by 3%, 16% respectively, for Fitlup, Fit2p)p, and by 5% 
for Fitp.p. Except the second case, which is ruled out 
by K~ data, the deviation from elastic FSI 

is mild. This reflects the fact that the color-suppressed 
rate is still small compared to Zl(*)+ 7 r“. The 
strong phases |i 5 i /2 — ^ 3 / 2 ! in Table p agree rather well 
with the directly extracted ones (Eg. (^)) as well as the 
elastic ones, and the validity of Eq. ( |^ as good estimates 
is born out. 

A principle motivation and interest in understanding 
color-suppressed B £)(*)o/j(*)o modes is its possible im¬ 
plication for B —!■ ATtt and tttt final states. These modes 
have been one of the focal points in B physics in recent 
years, because it provides rich probes of CP violation |2|] 
and possibilities for new physics. We note that the 
effect of 02 is rather subdued in these processes, but our 
picture of rescattering may still be realized, hence these 
processes provide more fertile testing ground for FSI. Ef¬ 
fects of FSI rescattering on ATtt, tttt final states have been 
discussed in the literature [|[ 011 . In particular, 

it has been stressed that large rescattering phases in 
Ktt Ktt and tttt —> tttt could have dramatic impact 
on such charmless final states. If phases analogous to 
6 ' ~ 50° can be realized, the and 7 r° 7 r° modes could 
get enhanced while tt+tt” mode suppressed. Direct CP 
rate asymmetries could soon be observed in Ktt modes, 
in particular in “pure penguin” K^tt~ mode, while for 
7 r° 7 r° and tt+tt” modes they could even reach 50-60%. 
The formalism is a straightforward extension from the 
one presented here. Rather than DH —> DH rescatter¬ 
ing, one now needs to study HH —s- HH rescattering. It 
is interesting to note that the factorization “sources” for 
all four Kn modes are sizable, unlike our present case 
where the mode is the singly large source. Thus, 

the cross-feed between channels would be different. Fur¬ 
thermore, the CP violating rate asymmetries provide ad¬ 
ditional leverage to check for the presence of FSI phases. 
In this sense the Ktt, tttt {PP in broader sense) system 
is richer than our present D^*^P case. It is interesting 
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that the physical picture of is still applicable with an 
additional annihilation rescattering term, due to possible 
final states consist of PP. Our study is underway and 
would be reported elsewhere. 

Finally, we mention a curiosity. Fit2£)p contains large 
“annihilation” rescattering, which runs against 
DfK~ data hence is ruled out. Fitlpp and Fitp.p, as 
well as Solution 1 of DV case, all had exchange rescat¬ 
tering far dominant over annihilation. Thus, Solution 2 
of DV case is the only one where the latter is sizable 
and dominant. It therefore has the distinct feature that 
is almost twice as large as D^u>, with Df K~ not 
much smaller. However, if this were realized, then one 
would expect —s- to be at factorization rate 

and much weaker than and D*^K°, which would 

be rather peculiar. We would therefore not be surprised 
if Solution 2 gets ruled out soon, and one might then 
conclude that rescattering is largely in terms of the clas¬ 
sic “charge exchange” type. This may also explain why 
KK modes are so far unseen. In this vein we wish to 
remark also that we have not exhausted the predictive¬ 
ness of our approach. For instance, one could generate 
“wrong charge” B ^ D K decays via D~Tr^ D K~ 

and D~Tr'^ D rescattering from Vub suppressed 
B ZIj 7 r°, decays. The rescattering matrix can 

be adapted from results presented here, but the latter 
decays have yet to be observed. 

In conclusion, we advocate in this work the possibility 
that the recently observed host of B ^ D^hP modes may 
be hinting at final state rescattering. In contrast to a sug¬ 


gestion of a larger and complex 02 , we extend the elastic 
FSI picture to quasi-elastic D^Pp 
D^Pp and DV DV rescattering, where P is the pseu¬ 
doscalar SU(3) octet, and V is the vector U(3) nonet. 
In this way we are able to accommodate D^*Pri, D^uj 
modes in a unified setting. For D^Pp modes, we find 
that data give rise to two rescattering phases 6 ' ~ 50° and 
6 ~ 20 °, where the need for a large phase comes about 
because of the strength of both the D^*pTr^ and D^*Prj 
modes. For DV modes, nonet symmetry reduces the 
number of physical phases to one, of order 20-30°. The 
emerging pattern is that of “charge exchange” rescatter¬ 
ing, rather than “quark annihilation”. We predict rather 
small B^ Di*'^~^K~, DfK^*'>~ and D^cj) rates, and 
£) 0 ^o ^ although in one solution one could have 

HOpO ~ 2 X D^uj and W D+K<^P-_^ 3 x 10^, 
which can be easily checked. We expect B^ 
and D^K*^ to be at 0.55 x 10“"^ and 0.2 x 10“"^, respec¬ 
tively, which is sizable. While these predictions can be 
tested experimentally, the Ktt, tttt charmless final states 
are even more promising, if FSI phases are as large as 
50°, because one then expects rather sizable direct CP 
asymmetries with a distinct pattern. 
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APPENDIX A: SU(3) DECOMPOSITION OF THE RESCATTERING MATRIX 


It is well known that Df and tt, K, 77 transform respectively as 3 and 8 under SU(3), 


/4 + 


D{3) = {D° D+ Dt ) , n(8) = 


V 


%/2 V 6 
7r“ 

K- 


K+ \ 


V 2 V6 


(Al) 




The 11(3) 0 n(8) can be reduced into a 3, a 6 and a 15, i.e. (see for example, |13[ ) 

r(3), = Din\, 


1 ^ ^ 7-1 TT^ 


r(15)}fe = DkB) + D,Bl - -5lDi^ - 


(A2) 


The SU(3) symmetry of strong interaction enforces the twenty-four by twenty-four scattering matrix having the 
following form 


15 6 3 

51/2 ^ ^ |r(T5); a) (T(T5); a| + ^ |r( 6 ); b) (r( 6 ); b\ + ^ |r( 3 ); c) (T(3); c|; (A3) 

a—1 b—1 c—1 

where |T(15);a), |T( 6 ); 6 ) and |r(3);c) are orthonormal SU(3) basis for the irreducible representations shown in the 
above equation. 
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Although the twenty-four by twenty-four rescattering matrix is diagonal in these basis, we may not need all of them 
in a realistic situiation. For a final state with given strangness and isospin (or electric charge) it can only rescatter 
to other final states having same quantum numbers. For later purpose, we giv e the explicit forms of these basis. By 
proper linear combining states within the same multiple, as shown in Eq. ( |A2f ), |T'(15);a) can be classified according 
to strangthness {S) and isospin (I) 


{S = 1,1=1) 

(^ = l,/=i) 

(5=1, J = 0) 
(^ = 0,/=^) 
(5 = 0,/=i) 

(5 = -l, J=l) 


\D+K+) + \D+Tr+) \D+K°)-\D°K+) - V2 \D+tt°) \D°K°) + \D+tt-) 

71 ’ 2 ’ 71 ’ 

\DtK+), iDtK^, 


\D+K°) + \D°K+) + V6\Dtm) 

272 


\D+tt+), l\D\+) + ^\D+n°), -^\D+7r-) 

2\D°tt+) - V2\D+Tr°) + 3V6\D+r)s) - 6\D+K°) 

iTl ’ 

2\D+tt-) + V2\D°Tr°) + 376111 %) - 6\D+K-) 

471 ’ 


\D+K^), 


\D+K-) + \D°K°) 

71 


111%-). 



Similarly |r(6),7 are 


(5=1, 7=1) : 


(*5 = 0,/=i) 


(5 = -1,7 = 0) : 


171+77+) - |71+7r+) iTl+TfO) - |71%+) -f 72|71+7r°) 

71 ’ 2 

2|71°7r+) - 72|71+7r0) - 76|71+7?8) - 2|71+770) 

4 ’ 

2|71+7r-) -f 72|71°7r0) _ 76|71%) - 2|71+77-) 

4 

|7i+a:-) - 1710770 ) 

71 ■ 


|71%°) - ITI+tt-) 

72 


The |r(3),c) are 


(5 = 0, 7=1/2) 


(5=1, 7 = 0) 


61710%) _ 3y2|71+7rO) -f 76|71+?78) -f 6|71+770) 

iTl 

6|71+7r-) -H 372|71°7r°) -f 76| jl°??8) + 6|il+77-) 

476 

3|71+77°) -f 3|71°77+) - 76|71+7?8) 

276 


(A4) 

(AS) 

(A6) 

(A7) 

(AS) 

(A9) 

(AlO) 

(All) 

(A12) 

(A13) 


(A14) 

(A15) 


With these basis, it is then straightfarward to obtain 5 ^A by using Eq. (|A3|) . As noted before we only need final 
states with same quantum numbers for rescattering. For example, 7 = 3/2 states can only appear in 15, hence we 
identify 153/2 as Jyg-. For 5 = — 1 and Q = 0, we only have and D'^K~ for rescattering. By using two neutral 

and 5 = — 1 states in 15 and 6, respectively, we immediately obtain 5 ^A ju Eq. (|7l| ). Similarly, by using 5 = Q = 0, 
7 = 3/2, 1/2 states in 15 and 7 = 1/2 states in 6 and 3, we obtain ju Eq. ([40D readily. 

We now turn to the DV case. The SU(3) decomposition of DV final states can be obtained by replacing tt, K and 
7/8 in the DP case by p, K* and a;8, respectively. However, there is another 3 as |T(3');7) = \D{3)iUJi). This 3' can 
mix with the previous 3 with a two by two symmetric (due to time reversal invarint) unitary matrix . We have 


57%e%^|r(15);a)(T(15);a| 

0=1 


e*^«^|r(6)7)(r(6)7H- ^ ^ |r(m);c)Z777(r(n);c|, 

^—1 m,n—3,3' c—1 


(A16) 
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FIG. 5: Geometric representation of B —> DP rescattering of Eq. (Bl). 

for DV rescatering matrix. Note that the symmetric mixing matrix can be parametrized by two phases and one 
mixing angle. We need four paramters, including three phase differences and one mixing angle, to discribe the DV 
FSI case. On the other hand by nonet symmetry we reduce them to only one parameter (in addition to an overall 
phase) but with two distinct solutions as discussed in Section IV. 


APPENDIX B: GEOMETRIC REPRESENTATION 


We give the geometric (triangular) representation of our results in this Appendix. For simplicity of presentation, 
we consider the leading FSI contribution and drop the 02 contribution. The triangular relation for the Dtt system, 
Eqs. (^) and ([^), and the FSI formula for Df K~, Eqs. (^) and (^, will then be reduced to 


-^3/2 + ■v/2Ai/2, 

= —2A3/2 + '\/ 2 Ai / 2 , 
= 3 A 3 / 2 , 

A3/2 = 


V2A,f2 = 4(1 + + 9e*®)(V2e*'*^/M(,2), 


16' 


V2AD0rj = ^(3 - 6e*^' + 3e*®)(V2e*'*3/M(/2)> 


1/2/ 


A 


DJK- 

V3 


= 4(-l - 2e*^' + 3e*®)(y2e*^3/^A{/2), 


'^■^ 1/2 ~ ^^3/2’ 


(Bl) 


= —+ '^^ 1/2 ~ *1 when 02 = 0. Eq. (^) still holds, and 


where the last equation follows from 

IA 3 / 2 I = 1 ^ 3^211 have IA 1 / 2 I < |A(/ 2 I 1° quasi-elastic rescattering. 

We illustrate the amplitudes of Fitlup, i.e. 5' = 47.8° and 9 = 24.8°, in Fig. |[ The D*P case is similar. We have 
chosen the cc-axis to coincide with A 3 / 2 . Since <5' and 6 are of the same sign, when considering only there is 

no need for a large angle, and 6 ' ^ 9 ^ 6 = S 1/2 — S 3/2 ^ 30° would have been good enough. However, because 9 has 
same sign as 6 ', as we can see from Fig. ||, the third phaser contributing to A^o^, i.e. 3e*®/16(e*^3/2.y2A:[^2)j turns 

back and tends to reduce We therefore need to start with a larger <5' to compensate. Thus, the relatively 

large phase 5' ^ 50° is driven by the strength of the measured B{D'^ri) and Note that in this case the S' 

and 9 phases compensate strongly for each other and A^+^_ is small compared to Apo^ and Apo^o. 

For Fit2 of DP case, 5'{= 17.1°) and 9{= —52.7°) are opposite in sign. This favors the generation of D'^rj since the 
effects of S' and 9 add to each other. However, we would need a large 9 phase to overcome the effect of S' to generate 
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Zl° 7 r° mode i.e. to account for |i5| = \5i/2 — ^ 3 / 2 ] ~ 30°. Otherwise, we will have too small a A^Ot^o amplitude. Same 
as the previous case, a large phase of order 50° is needed to account for the strength of both the Z)° 7 r° and D^rj 
modes. Inspecting the generation of however, we see that the effects of 5' and 6 now add to each other, and 

would generate too large a B{DfK~) that is already ruled out by data, hence the case is not plotted. 

We also refrain from plotting the case for DV modes for the following reasons. First of all, because of relatively 
weaker rescattering, the 02 effect is more prominent than in D^*'>P case. Second, some discussion is already given in 
Sec. IV.B, where comparison is made between FSI and 02 contributions. Decomposing into A 1/2 and A ^/2 amplitudes 
does not make the case clearer. ,_ 
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